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The LHC (CMS) discovery potential of the e^fi^ + Ej? tss signature in the search 
for supersymmetry and lepton flavour violation in neutralino decays is studied. A de- 



< 

in 



tailed study is done for the CMS test points LM1-LM9. It is shown that for the point 
LM1 it is possible to detect lepton flavour violation in neutralino decays with lepton 



flavour violating branching Br{\2 ~^ ^ ±e ^Xi) — 0.04Sr(x2 — ► e+e XiiM^A 4 Xi) f° r 
I — . an integral luminosity 10/6 . A discovery potential in the mSUGRA-SUSY scenario 



One of the goals of the Large Hadron Collider (LHC) 1] is the discovery of supersymmetry 



with tan/? = 10, sign(fi) = + in the (mo, m^j) plane using the e^p^ + E™ 
' signature is determined. 

a: 

I. INTRODUCTION 

> 
• r- 1 

X 

(SUSY). The squark and gluino decays produce missing transverse energy from lightest stable 
superparticle (LSP) plus multiple jets and isolated leptons [lj. One of the most interesting 
and widely discussed signatures for SUSY discovery at the LHC is the signature with two 
opposite charge and the same flavour leptons [2|: + E™ lss . The main reason of such 

interest is that neutralino decays into leptons and LSP \\ ~ * l + l~Xi contribute to this 
signature and the distribution of the invariant mass m inv (l + l~) has the edge structure 
3] that allows to determine some combination of the SUSY masses. 



*Electronic address: Iouri.Andreev@cern.ch 
^Electronic address: Serguei.Bitioukov@cern.ch 
^Electronic address: Nikolai. Krasnikov@cern.ch 
§ Electronic address: Alexander. Toropin@cern.ch 



2 



The signature e ± /i T + E™ lss , can also be produced when the decays into r pair which 
is only relevant at large tan/3 JJ|. Also as it is shown in Ref. 3[ at the level of CMS JET 
simulation the use of e ± /i =F + E™ ss signature for large tan j3 allows to obtain nontrivial 
information on parameters of the decay X2 ~ > — > ttx\ — > e yL^y^yvvv. On the other 
hand, the e^y^ + E™ lss (with an arbitrary number of jets), can be used for the detection 
of lepton flavour violation in slepton decays 

B-Q atthe LHC. 

In the Minimal Supersymmetric Model (MSSM) supersymmetry is broken at some 
high scale M by generic soft terms, so in general all soft SUSY breaking terms are arbitrary 
which complicates the analysis an d sp oils the predictive power of the theory. In the Minimal 
Supergravity Model (mSUGRA) the universality of the different soft parameters at the 
Grand Unified Theory (GUT) scale Mqut ~ 2- 10 16 GeV is postulated. Namely, all the spin 
zero particle masses (squarks, sleptons, higgses) are postulated to be equal to the universal 
value m at the GUT scale. All gaugino particle masses are postulated to be equal to the 
universal value mi 12 at GUT scale. Also the coefficients in front of quadratic and cubic 
SUSY soft breaking terms are postulated to be equal. The renormalization group equations 
are used to relate GUT and electroweak scales. The equations for the determination of 
a nontrivial minimum of the electroweak potential are used to decrease the number of the 
unknown parameters by 2. So the mSUGRA model depends on five unknown parameters. At 
present, the more or less standard choice of free parameters in the mSUGRA model includes 

masses depend on these parameters. 

The goal of this work is to search for the possibility to detect SUSY and lepton flavour 
violation (LFV) using the e^p^ + E™ lss signature at the LHC for the Compact Muon 
Solenoid (CMS) detector at the level of full detector simulation. For specific calculations 
the mSUGRA model is used. 

The organization of the paper is the following. Section |H] describes some useful technical 
details of performed simulations. In Section ITTTl the backgrounds and cuts used to suppress 
the backgrounds are discussed. Section HXl contains the results of numerical calculations 
concerning the possibility to detect SUSY using the e ± fi T + E™ lss signature. In Section IVl 
the prospects of the detection lepton flavour violation in the neutralino decays is studied. In 
Section IVll the influence of the systematic uncertainties on the value of the signal significance 
is discussed. Section I VI II contains concluding remarks. 
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II. SIMULATION DETAILS 

The coupling constants and cross sections in the .ead.ng order (LO) applanation for 
SUSY processes and backgrounds were calculated with ISASUGRA 7.69 |l6j, PYTHIA 



6.227 |12( and CompHEP 4.2pl For the calculation of the next-to-leading order (NLO) 
corrections to the SUSY cross sections the PROSPINO |3] code was used. For considered 
signal events and backgrounds the NLO corrections are known and the values of NLO cross 
sections (or /c-factors) were used for normalization of the numerical results. 

Official datasets (DST) production was used for the study of CMS SUSY test points 



0, |2l| (LM1, LM4, LM5, LM9), lepton flavour violation for the point LM1 and of back- 
grounds (tt, ZZ, WW, Wt, Zbb, DY2r). The ISASUGRA 7.69 + PYTHIA 6.225 codes 
were used in official production. The full detector simulation was made with OSCAR_2_4_5 
or OSCAR_3_6_0 Isj codes. Digitization was made with ORCA_7_6_l, ORCA_8_5_0 or 



ORCA_8_7_l |22j codes. 

For other CMS SUSY test points (LM2, LM3, LM6, LM7, LM8) and WZ background 
;he events were generated with ISASUGRA 7.69 + PYTHIA 6.227 codes and CMKIN_4_3_1 
22] ] was used as an interface program. The detector simulation and hits production for 
the test points (LM2, LM3, LM6) and WZ background were made with OSCAR_3_6_5 
and for digitization ORCA_8_7_3 was used. To study the two test points LM7 and LM8, 
background Z+jet and to prepare the CMS discovery plot, the CMS fast simulation program 
FAMOSAAD |22| was used. 

The pile-up for the signal events are not taken into account, but backgrounds in DSTs 
were produced with pile-up corresponding to 2 x 10 33 cm _2 s _1 luminosity. 

The reconstructed electrons and muons were passed through packages defining lepton 
isolation criteria. For each electron and muon the following parameters are defined: 

• Trackl solation is a number of additional tracks with > 2 GeV/c inside a cone with 



R = ^Ar] 2 + A$ 2 < 0.3 around the lepton. 

Calol solation is a ratio of energy deposited in the calorimeters (electromagnetic 
(ECAL) + hadronic (HCAL)) inside a cone with R = 0.13 around given track to 
the energy deposited inside a cone with R = 0.3. 

HEratio is defined as a ratio of energy deposited in the HCAL inside a cone with 
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R = 0.13 to the energy deposited in the ECAL inside the same cone. 

• EPratio is a ratio of energy deposited in the ECAL inside a cone with R = 0.13 to 
the momentum of the reconstructed track. 

Data with reconstructed electrons, muons, jets and missing energy were stored into ROOT 
(23! files for the final analysis. 



The official datasets used for these analysis were processed with CRAB 
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III. SIGNAL SELECTION AND BACKGROUNDS 

The SUSY production pp — > qq , gg, qg with subsequent decays 

<?^<?'xf 2 > (1) 
9 -> Qq'xffl, (2) 
Xt,2 -> X?e + (//>, (3) 
xr, 2 XiM~(0^> (4) 
leads to the event topology e ± /i =F + E™ lss . Note that in the MSSM with lepton flavour 
conservation neutralino decays into leptons X234 — > l + l~Xi = e 'A i ) do not contribute into 
this signature and contribute only to the + E™ lss signature. The main backgrounds 

which contribute to the e ± ^ events are: tt, WW, WZ, ZZ, Wt, Zbb, DY2r and Z+jet. It 
is found that tt is the largest background and it gives more than 50% contribution to the 
total background. The following NLO values for the main background cross sections [3, 13 
are used (Table HJ). 

In the final analysis the events with the following isolation criteria for electrons were used: 
Tracklsolation < 1.0, C alol solation > 0.85, 0.85 < EPratio < 2.0, HEratio < 0.25. The 
same criteria for muons were the following: Trackl solation < 1.0, C alol solation > 0.50, 
EPratio < 0.20, HEratio > 0.70. These numbers were adjusted by studying electron and 
muon tracks in the process pp — > WW — ► 21. 

The selection cuts are the following: 

• cut on leptons: p l ^ pt > p l ^, pt, °, \n\ < 2.4, lepton isolation within AR < 0.3 cone 

• cut on missine; transverse energy: Ej? tss > Erp ' . 
Where p l ^ pt '° and E™ %ss '° are corresponding thresholds. 
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TABLE I: The main background cross sections (in pb). 



Process 


(7 T D 


(jATTCi 
^ 1 V LjKJ 


tt 


505 


830 


WW 


70 


117 


wz 


27 


50 


zz 


11 


16 


Wt 


30 


62 


zbb 


790 


1580 


Z+jet, ckin(3) 


= 100 240 


274 


DY2r 


39600 





A. Trigger selection 

The events are required to pass the Global Level 1 Trigger (LI) 27\ and the High Level 
Trigger (HLT) [28]. The events have to pass at least one of the following triggers: single 
electron, double electron, single muon, or double muon. The used cut on leptons is more 
stringent than the cuts used in the HLT for these triggers. 

IV. USE OF THE e ± u T + EW iss SIGNATURE FOR THE SUSY DETECTION 



20 



are chosen 



The possibility to detect SUSY using the CMS test points LM1 - LM9 
for the detailed study of SUSY detection at CMS is investigated in this section. This study 
is based on the counting the expected number of events for both the SM and the mSUGRA 
models. The parameters of the CMS test points LM1 - LM9 are given in Table ITT1 

For the point LM1 (the point LM1 coincides with the post-WMAP point B Q) the 
distributions on p l ^ pt , EJp tss and m inv (e + fi~ + e~/i + ) for both background and signal events 
are shown in Figs HE] 

It was found that the set of cuts with p^ pt > 20 GeV/c, EJp iss > 300 GeV is close to the 
optimal set (the highest significance with the best signal/background ratio). The results for 
the luminosity £ = 10 fb~ are presented in Table ILTTl 

For other CMS SUSY test points LM2 - LM9 the results with the same set of cuts 
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TABLE II: The parameters of the CMS test points. 



Point 


m (GeV) 


m i/2(GeV) tan/3 sign(/i) 


A 


LM1 


60 




i n 

1U 


i 

T 





LM2 


185 


OOU 


DO 


1 
I 





LM3 


330 




on 


1 
1 





LM4 


210 


285 


10 


+ 





LM5 


230 


360 


10 


+ 





LM6 


85 


400 


10 


+ 





LM7 


3000 


230 


10 


+ 





LM8 


500 


300 


10 


+ 


-300 


LM9 


1450 


175 


50 


+ 






are presented in Table IIVI The significances definitions in this table are the followin 

S c i2 = 2(VN S + N B - v^Vfl) Q and S cL = ^2((N S + N B ) ln(l + |f) - N s ) El 



TABLE III: The expected number of events for backgrounds and for signal at the point LM1, 
C = 10 fb- 1 , e ± f^ + EJp iss signature. 



Process 


2 isolated leptons, p l ^ pt > 20 GeV/c E^ iss 


> 300 GeV 


tt 


39679 


79 


WW 


4356 


4 


wz 


334 


2 


TIL 


38 





Wt 


3823 


2 


Zbb 


315 





Z+jet 


1082 


6 


DY2r 


7564 





SM backj 


ground 57191 


93 



LM1 Signal 



1054 



329 
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TABLE IV: The number of signal events and significances for cut set with p^ p > 20 GeV/c and 
EJp iss > 300 GeV for £ = 10 fb' 1 , signature e ± /i =F + EJp iss . The number of the SM background 
events N B = 93 (see Table EJ . 



Point 


TV events 




c 

<->cL 


T A T1 


329 




z4.y 


T A TO 


94 


0.1 


0.0 


t ivn 

IjlvlO 


AD9 


9^ 9 


9Q 9 


LM4 


301 


20.4 


23.1 


LM5 


91 


7.8 


8.3 


LM6 


222 


16.2 


18.0 


LM7 


14 


1.4 


1.4 


LM8 


234 


16.9 


18.8 


LM9 


137 


11.0 


11.9 



It was found from the Tables ITLTllIVI that for the point LM1 the significances are S c \2/S c l = 
21.8/24.9 for the e ± /i T + E^ iss signature. 

The supersymmetry discovery potential for the mSUGRA model with tan/5 = 10, 
signal) = + in the (mo, rnx/2) plane (generalization of the point LM1) using the CMS 
fast simulation program FAMOS AAD 22j was also studied. The CMS discovery potential 



contours for C = 1, 10 and 30 fb 1 for the signature e ± /i =F + E™ lss are shown in FiglHl 



A. The comparison of the FAMOS and the full simulation for the point LM1 

The results obtained with the FAMOS code were compared with the full simulation 
results. The test point LM1 was used as the comparison object. The distributions on p l ^ pt 
and E?p iss of the LM1 signal for both FAMOS and ORCA are shown in FigsEEl 

From the Figs. 15161 it is possible to conclude that the full and the fast simulation distri- 
butions on p l ^ pt and E™ zss for the point LM1 are in reasonable agreement. 
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V. SEARCH FOR LEPTON FLAVOUR VIOLATION IN THE NEUTRALINO 

DECAYS 

In the MSSM the off-diagonal components of the slepton mass terms violate lepton flavour 
conservation. As it has been shown in Refs. it is possible to look for lepton flavour 
violation at supercolliders through the production and decays of the sleptons. For the LFV 
at the LHC, one of the most promising processes is the LFV decay of the second neutralino 
|il [ill xt M ~~* Xi tt' -i where the non zero off-diagonal component of the slepton mass 
matrix leads to different lepton flavours in the final state. This mode is more sensitive to 
LFV compared to the direct Drell-Yan production of sleptons since the second neutralino x*2 
can be copiously produced through the cascade decays of squarks and gluinos Q] . By using 
the above mode, LFV in e — jl mixing has been investigated in Refs.j^.fiol] at a parton 
model level for a toy detector. In this section the perspectives of LFV detection in CMS on 
the base of full simulation of both signal and background is studied. To be specific, the test 
point LM1 is studied. The signal of LFV Xz decay is the two opposite-sign leptons (e + /i~ 
or e~fi + ) in the final state with the characteristic edge structure. In the limit of lepton 
flavour conservation, the process ~ * H ~ > HXi h &s the edge structure for the distribution 
of lepton-pair invariant mass m u and the edge mass m™" a: (7 + Z _ ) is expressed by slepton mass 
mr and neutralino masses m^o as follows |3|: 

777? m 2 ~ 

(mz x (i + nr = ™|o(i - -f )(i - -*o (5) 

The SUSY background for LFV comes from uncorrelated leptons from different squark 
or gluino decay chains. The SM background comes mainly from 

it -> bWbW -> blbivv. (6) 

It should be stressed that the signature with e^y^ in the absence of LFV do not have the 
edge structure for the distribution on invariant mass m m „(e ± 1 u T ). As the result of LFV the 
edge structure for e^y^ events arises too. Therefore the signature of LFV is the existence 
of an edge structure for the e ± fi T distribution. The rate for a flavour violating decay is 
determined by nonzero k, where 
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In this paper the observability of LFV for the point LM 1 was studied. For this purpose 
a special sample of events with 100% LFV was prepared. For k = 0.10, k = 0.25 the 
distributions of the number of e ± fi T events on invariant mass rrii nv (e (FigEJ) clearly 
exhibit the edge structure, i.e. the existence of the lepton flavour violation in neutralino 
decays. 

It appears that for the point LM1 the use of an additional cut 

m ira (e ± /i T ) < 85 GeV (8) 

reduces both SM and SUSY backgrounds and increases the discovery potential in the LFV 
search. In this case the cuts p l T Pt > 40 GeV/c and E™ ss > 200 GeV give the best result. 
For instance, it was found that N S Mb g = 19 and N S usYb g = N LM i signa i = 71 for the e ± fi T 
signature (see Table ILTTl for comparison where NsMbg = 93 and NiMisignai = 329). 

It was found that for the point LM1 in the assumption of exact knowledge of backgrounds 
(both SM and SUSY backgrounds) for the integral luminosity of £ = 10/6 -1 it would be 
possible to detect LFV at a 5a level in \2 decays for k > 0.04. 



VI. INFLUENCE OF THE SYSTEMATIC UNCERTAINTIES ON THE 

SIGNAL SIGNIFICANCE 

In this Section the influence of systematic uncertainties on the value of signal signifi- 
cance in the case of SUSY detection is estimated. The systematic uncertainties in the signal 
significance calculation include the experimental selection uncertainty of the background 
events, luminosity uncertainty, and the theoretically calculated uncertainties of the tt, WW 
and other backgrounds. Theoretically calculated uncertainty in background cross sections 
consists of the uncertainty related with inexact knowledge of parton distribution functions 
(PDF) and higher order corrections to the NLO cross sections 3^]. There are several exper- 
imental uncertainties related with lepton identification, missing energy and luminosity. In 
accordance with Ref . the systematic error related with the lepton identification is 3%, 
the systematic error related with the missing energy is 2% j^J. 

It was found that the experimental uncertainties in the number of background events 
Nb related with the missing energy and the lepton identification lead to 10% and 0.5% 
uncertainties in the number of the Nb respectively. 
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The systematic uncertainty in the luminosity is 5% 35]. The total 5% uncertainty in 
luminosity leads to 5% uncertainty in the number of background events. 

In the studying signature the tt background dominates. The PDF uncertainty of tt cross 
section is equal to 5% and the uncertainty due to unknown higher order corrections to 
the NLO background cross section is equal to 10%. In the assumption that the systematic 
uncertainties are added quadratically it was found that the overall uncertainty in the number 
of background events is about 16%. 

Following the prescriptions of the CMS PRS group the influence of the systematic uncer- 
tainties on signal significance using the program for calculations of significance S_cP from 



Ref. 



was calculated. The results are shown in Table IV! 



TABLE V: The dependence of signal significance S c \2 on background uncertainty for the used set 
of cuts and C = 10 for signature e ± ^z T + EJp iss . 



Point 


0% 


10% 


20% 


LM1 


21.8 


15.7 


10.0 


LM2 


8.1 


5.8 


3.7 


LM3 


25.2 


18.1 


11.6 


LM4 


20.4 


14.7 


9.4 


LM5 


7.8 


5.6 


3.6 


LM6 


16.2 


11.7 


7.5 


LM7 


1.4 


1.0 


0.6 


LM8 


16.9 


12.1 


7.8 


LM9 


11.0 


7.9 


5.1 



The discovery potential decreases with the increase of the background uncertainty but 
nevertheless the results are rather robust. In particular, it was found that for 0, 10, 20% 
background uncertainties it is possible to detect at the point LM1 lepton flavour violation 
in neutralino decays for k > 0.04, 0.043, 0.051 respectively. 
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VII. CONCLUSION 

In this paper the possibility to detect SUSY and lepton flavour violation at the LHC 
(CMS) using the signature e ± /i =F + E™ lss was studied. This signature allows to discover 
both SUSY and LFV in the neutralino decays. It was found that for the CMS test points 
LM1 - LM9 for integral luminosity C = 10 fb~ x it is possible to discover SUSY for all LM 
points except for LM7. For tan/5 = 10, A = 0, sign(fi) = + the discovery contours for 
C = 1, 10 and 30 fb~ l in the (m , m^) plane were determined. The possibility to look for 
lepton flavour violation in neutralino decays was also studied. It was found that for the point 
LM1 it is possible to detect lepton flavour violation provided the lepton flavour violating 
branching Br(x2 ~~ * A t±e=F Xi) > 0.04£>r(x2 — > e + e~Xi, f° r £ = 10 fb" 1 . The 

results for the significances are rather robust under the inclusion of reasonable systematic 
uncertainties. 

It should be stressed that the signature e ± fi T + E™ lss is less "powerful" from the point 
of view of SUSY discovery that say the signature (n > 3 jets) + E™ lss 0, however it is 
important to detect "new physics" using simultaneously different signatures that increases 
the credibility of the results. 
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FIG. 1: The p l ^ pt distribution after selection of two isolated leptons with p l ^ pt > 20 GeV/c. The 
both leptons are plotted. 
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FIG. 2: The E™ lss distribution after selection of two isolated leptons with p^ pt > 20 GeV/c. 
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FIG. 3: The p^ p distribution after selection of two isolated leptons with p^ p > 20 GeV/c and 
EJp iss > 300 GeV. 
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FIG. 4: The invariant two lepton mass distribution after selection of two isolated leptons with 
p l * pt > 20 GeV/c and EJp iss > 300 GeV. 
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FIG. 5: The p l ^ pt distribution after selection of two isolated leptons with p l ^ pt > 20 GeV/c for 
ORCA and FAMOS. The plots are normalized to the numbers of dileptons. 
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FIG. 6: The E™ tss distribution after selection of two isolated leptons with p£ pt > 20 GeV/c for 
ORCA and FAMOS. The plots are normalized to the numbers of dileptons. 



20 




FIG. 7: The distribution of two lepton invariant mass after selection of two isolated different-flavour 
leptons with p l ^ pt > 20 GeV/c and > 300 GeV for flavour violation parameter k = 0.10 and 

k = 0.25. Superimposed curves are fits to the invariant mass distribution for 100% LFV. 
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FIG. 8: e + fi + e /i + discovery plot for tan/? = 10, sign(fj,) = +, A = 0. Selected two isolated 
leptons with p l ^ pt > 20 GeV/c and Ef iss > 300 GeV. Calculations are made for S cl2 . 



